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Preface 


This  thesis  is  a  closer  look  at  an  atmospheric  radio  noise 
(ARN)  model  proposed  by  Steven  Hettinger,  a  former  AFIT  student. 
Hettinger  found  that  the  existing  models  of  ARN  are  not  tractable 
through  existing  receiver  structure  and  that  higher  order  statistics 
are  not  available  for  these  models.  The  proposed  model  is  simple 
enough  to  be  tractable  and  yet  complex  enough  to  model  many  types  of 
ARN.  In  this  report,  the  statistics  of  the  model  are  derived,  an 
optimal  receiver  for  communication  through  this  noise  is  computed,  and 
the  accuracy  of  the  model  is  examined  by  comparing  the  first  order 
statistics  of  the  model  to  measured  first  order  statistics. 

I  would  like  to  thank  my  advisor.  Major  Joseph  Carl  for  his 
support  and  guidance  throughout  this  research.  I  would  also  like  to 
thank  Captain  Stanley  Robinson  for  his  advice  in  many  areas  and 
Captain  Pedro  Rustan  for  his  aid  in  understanding  the  nature  of 
lightning.  I  must  also  give  special  thanks  to  my  wife  Rosann  for  her 
support  throughout  the  project. 

John  F.  Stach 
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Abstract 

An  investigation  of  low  frequency  atmospheric  radio  noise 
(ARN)  indicates  that  the  return  stroke  from  lightning  discharges  is 
the  major  source  of  the  noise.  A  model  of  ARN  is  then  presented  which 
is  based  on  the  return  stroke. 

The  model  consists  of  the  sum  of  two  marked  Poisson  processes. 
The  sum  of  marked  Poisson  processes  can  be  thought  of  as  a  single 
marked  Poisson  process  with  a  transformed  rate  and  a  mixture  density 
on  the  marks.  The  complete  statistical  description  of  this  process  is 
then  derived  and  a  minimum  probability  of  error  processor  is  designed 
using  the  ARN  model  as  the  noise. 

The  theoretical  amplitude  probability  distribution  (APO)  is 
then  derived  for  the  output  of  a  quadrature  envelope  detector  using 
the  ARN  model  as  an  input.  The  theoretical  APD  curves  were  then  com¬ 
pared  to  a  measured  APD  curve  from  CCIR  report  322  and  two  measured 
APD  curves  from  an  article  by  K.  Furutsu  and  T.  Ishida.  The  linear 
mean  square  error  (MSE)  between  the  theoretical  and  measured  curves 
when  plotted  on  Rayleigh  paper  is  .018  and  .015  for  the  Furutsu  and 
Ishida  data,  and  .018  for  the  CCIR  data. 

Based  on  the  first  order  statistics,  the  model  is  found  to  be 
a  good  representation  of  low  frequency  ARN.  It  is  therefore  recom¬ 
mended  that  the  model  be  used  to  evaluate  receiver  performance  in  the 
low  frequency  channel  and  that  a  more  thorough  validation  study  be 
done  including  a  validation  of  higher  order  statistics. 


VALIDATION  AND  RECEIVER  DESIGN 
FOR  A  RANDOM  POINT  PROCESS  MODEL  OF  ATMOSPHERIC 
RADIO  NOISE 

I.  Introduction 

The  increasing  demand  for  communications  has  created  an 
increasing  demand  for  usable  frequency  bands.  The  low  frequency 
spectrum  (D.C.  to  30  MHz)  is  usable,  but  it  has  a  type  of  noise  that 
is  different  than  in  the  high  frequency  spectrum.  In  the  low  fre¬ 
quency  channel,  atmospheric  noise  is  frequently  the  limiting  factor  on 
the  performance  of  communication  systems  (Ref  15).  It  is  necessary  to 
know  detailed  statistics  of  the  noise  in  order  to  calculate  the  chan¬ 
nel  capacity  or  error  rate  of  a  noisy  channel.  Also,  receiver  perfor¬ 
mance  can  be  analyzed  if  a  model  with  accurate  statistics  can  be 
traced  through  a  mathematical  representation  of  the  receiver. 

Finally,  an  accurate  model  for  which  the  complete  statistics  are  known 
will  allow  an  optimal  Bayes  receiver  to  be  designed  for  the  channel. 

In  this  report,  a  model  of  low  frequency  atmospheric  noise  (ARN)  is 
examined.  The  purposes  of  this  report  are  to  (1)  determine  the  vali¬ 
dity  of  the  model,  (2)  to  determine  the  possibility  of  using  the  model 
to  compute  the  performance  of  low  frequency  receivers  and  (3)  to 
determine  the  possibility  of  designing  optimal  receivers  for  the  low 
frequency  channel. 
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This  report  will  be  divided  into  four  major  sections.  (1) 
First,  the  causes  of  ARN  are  examined.  The  main  source  of  low  fre¬ 
quency  ARN  is  found  to  be  the  return  stroke  of  cloud  to  ground 
lightning  discharges.  (2)  Then,  a  sum  of  marked  Poisson  processes  is 
used  to  model  the  ARN.  It  is  found  that  the  resulting  process  can  be 
written  in  terms  of  a  single  Poisson  process  with  a  mixture  density  on 
the  marks.  The  minimum  probability  of  error  processor  is  then 
derived.  (3)  A  validation  study  is  then  done  to  determine  the 
accuracy  of  the  model.  Measured  amplitude  probability  distributions 
(APD)  of  the  output  of  an  envelope  detector  are  then  compared  to  the 
theoretical  APD  curves  of  the  model.  The  measured  APD  curves  are 
taken  from  CCIR  Report  322  (Ref  22)  and  an  article  written  by  Furutsu 
and  Ishida  from  Japan  (Ref  3).  The  mean  square  error  between  theore¬ 
tical  and  measured  APD  curves  plotted  on  Rayleigh  paper  is  found  to  be 
.018  and  .015  for  two  typical  curves  from  the  Furutsu  and  Ishida 
article,  and  .018  for  a  typical  curve  from  CCIR  Report  322.  (4) 
Finally,  some  conclusions  about  the  model  will  be  made  along  with  the 
recommendations  that  the  model  be  more  thoroughly  validated  and  that 
the  model  be  used  for  perfor-mance  analysis  of  receivers  for  the  low 
frequency  channel. 

The  next  section  describes  the  physical  processes  involved  in 
ARN  and  some  of  the  properties  of  ARN  that  will  affect  the  parameters 
of  an  ARN  model. 
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II.  Background 


There  are  several  types  of  noise  than  can  interfere  with  low 
frequency  communications.  In  the  low  frequency  spectrum,  man-made, 
galactic,  and  atmospheric  noise  are  the  typical  types.  Man-made  noise 
comes  from  power  lines,  machinery,  or  transmitters.  Propagation  of 
man-made  noise  other  than  transmitters  is  principally  through  power 
lines  and  by  ground  wave  (Ref  22:9).  Galactic  noise  power  figures  are 
available,  and  can  also  influence  low  frequency  noise  depending  on  the 
frequency  and  directional  properties  of  the  antenna  (Ref  22:11).  It 
is  important  to  mention  man-made  and  galactic  noise  because  of  the 
situational  factors  that  could  allow  these  types  of  noise  to  dominate 
a  received  waveform.  However,  the  articles  used  to  verify  the  model 
to  be  presented  attempt  to  isolate  the  atmospheric  noise  and  assume 
that  the  atmospheric  noise  predominates  at  the  frequencies  of  interest 
(Refs  3,  22).  Because  of  this,  the  model  is  designed  to  represent 
only  the  atmospheric  noise. 

At  any  given  time,  there  are  many  lightning  storms  around  the 
world.  The  total  electromagnetic  radiation  that  reaches  a  receiver 
from  these  storms  will  be  referred  to  as  atmospheric  noise.  Lightning 
is  the  main  source  of  atmospheric  noise.  Lightning  can  be  defined  as 
a  transient,  high  current  discharge.  Lightning  occurs  when  some 
region  of  the  atmosphere  attains  an  electric  charge  sufficiently  large 
that  the  electric  fields  due  to  the  charge  cause  electrical  breakdown 
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of  the  air.  Uman  (Ref  16)  treats  the  subject  of  lightning  thoroughly 
in  his  book  which  is  used  as  a  reference  in  the  following  discussion. 
Lightning  can  occur  within  a  cloud  (cloud  discharges),  between  two 
clouds  (cloud  to  cloud  discharges),  between  a  cloud  and  the  earth 
(ground  discharges),  or  between  a  cloud  and  the  surrounding  air  (air 
discharges).  The  type  most  people  are  familiar  with  are  ground  dis¬ 
charges.  A  total  ground  discharge  is  called  a  flash.  A  ground  flash 
is  composed  of  one  or  more  separate  strokes  in  the  same  or  separate 
ionized  paths.  Each  stroke  lasts  for  milliseconds  and  the  time  inter¬ 
val  between  strokes  is  roughly  50  msec.  The  flash  begins  when  the 
regions  between  some  of  the  positive  and  negatively  charged  areas  of  a 
cloud  achieve  a  potential  greater  than  the  ionization  potential.  A 
cloud  potential  stepped  leader  then  moves  toward  the  ground  in  roughly 
50  m  steps  with  a  pause  time  of  about  50  usee  between  each  step. 
Several  branches  may  propogate  toward  the  ground  until  one  or  more 
brings  cloud  potential  near  enough  to  the  ground  for  a  return  stroke 
to  occur.  When  the  stepped  leader  gets  near  the  ground,  the  resulting 
field  causes  a  column  of  ground  potential  charge  to  extend  from  the 
ground  to  the  leader  tip.  This  column  of  charge,  called  a  return 
stroke,  travels  with  a  velocity  from  one-tenth  to  one-third  the  speed 
of  light  up  the  ionized  path  that  the  leader  created,  reaching  the 
cloud  base  in  about  70  usee.  The  return  stroke  current  typically 
reaches  about  10  to  200  Kamps  in  a  few  microseconds  and  falls  to  half 
its  peak  value  in  about  20  to  60  usee.  Currents  on  the  order  of  hun¬ 
dreds  of  amperes  may  continue  to  flow  for  several  milliseconds.  If 
additional  charge  is  available  in  the  cloud  after  the  return  stroke. 
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the  flash  may  contain  additional  strokes.  These  strokes  usually  begin 
with  a  dart  leader  propagating  from  the  cloud  to  the  earth  along  the 
already  ionized  path  at  a  speed  of  about  2x10^  m/sec,  which  is  a 
magnitude  faster  than  the  stepped  leader.  The  dart  leader  is  then 
followed  by  another  return  stroke.  If  the  ionized  path  is  weakened 
after  the  return  stroke  by  wind  or  a  time  lapse  of  more  than  100  msec, 
then  a  dart-stepped  leader  may  occur  which  alternates  between  the 
stepped  and  dart  leader  as  it  moves  toward  the  earth. 

Cloud  discharges  are  the  most  frequently  occuring  form  of 
lightning.  They  generally  neutralize  about  the  same  amount  of  charge 
as  a  ground  discharge,  but  are  characterized  by  relatively  smooth 
field  changes  as  opposed  to  the  ground  return  stroke.  Also,  the  paths 
of  cloud  discharges  are  found  to  be  more  nearly  horizontal  than  verti¬ 
cal.  Although  these  differences  may  alter  their  statistics,  the  model 
of  ARN  will  not  distinguish  between  them  for  simplicity. 

Propagation 

Another  important  variable  in  the  study  of  atmospherics  is  the 
effect  of  propagating  medium  between  the  flashes  and  the  receiver. 

The  characteristics  of  the  medium  are  a  function  of  the  frequency  of 
the  emitting  source,  the  propagation  distance,  and  the  reflective  pro¬ 
perties  of  the  earth  and  ionosphere  (Ref  12).  The  two  common  modes  of 
propagation  of  low  frequency  radiation  are  the  groundwave  and  wave¬ 
guide  modes.  The  groundwave  has  little  effect  on  the  receiver 
measurements  because  the  antennas  used  to  collect  ARN  data  are  placed 
well  above  the  ground  to  avoid  groundwave  effects.  The  waveguide  mode 
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results  because  of  the  complex  conducting  properties  of  the  ionosphere 
and  the  earth.  The  waveguide  mode  of  propagation  allows  radiation 
from  distant  storms  to  reach  the  receiver.  The  properties  of  atmos¬ 
pheric  waveguide  propagation  have  been  summarized  by  Uman  (Ref  16), 
Harth  (Ref  5),  and  many  others. 

Power  Spectrum 

The  leader  and  return  stroke  affect  a  different  portion  of  the 
low  frequency  spectrum.  The  peak  amplitude  for  the  frequency  spectrum 
of  the  return  stroke  generally  occurs  between  5  and  10  KHz.  The  peak 
amplitude  of  the  frequency  spectrum  for  the  leader  stroke  is  typically 
between  3  and  30  KHz.  From  about  10  KHz  to  100  MHz,  the  peak  electric 
field  intensity  of  vertical  polarized  radiation  from  lightning  varies 
approximately  inversely  with  frequency  (Ref  16). 

In  summary,  the  characteristics  of  ARN  are  sensitive  to  the 
frequency  of  operation  and  the  location  of  the  receiver.  Since  the 
sun  affects  the  ionosphere,  the  characteristics  of  ARN  are  also 
affected  by  seasonal  changes  and  the  time  of  day.  To  represent  ARN 
with  a  model,  the  model  must  have  sufficient  parameters  to  account  for 
each  of  these  factors.  In  the  next  chapter,  a  model  will  be  presented 
which  is  based  on  the  physical  processes  discussed  in  this  chapter. 
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III.  Point  Process  Model 


An  accurate  noise  model  with  a  complete  statistical  descrip¬ 
tion  is  very  useful  in  receiver  design  and  evaluation  as  well  as  in 
other  aspects  of  electronic  communication.  Besides  providing  a  complete 
description,  the  model  must  also  be  simple  enough  so  that  it  can  be 
traced  through  a  mathematical  representation  of  a  receiver.  Such  a 
model  would  be  a  useful  tool  in  the  design  and  evaluation  of  receivers 
that  must  operate  in  that  channel.  In  addition,  the  model  must  be 
validated  to  insure  that  the  statistics  of  the  model  are  in  agreement 
with  the  statistics  of  the  actual  measured  noise  process.  In  this 
section,  a  model  is  proposed  for  which  the  complete  statistical  des¬ 
cription  is  derived  and  a  minimum  probability  of  error  receiver 
designed.  The  accuracy  of  the  model  will  be  investigated  in  the  next 
chapter. 

Development  of  the  Proposed  Model 

Atmospheric  noise  consists  of  a  highly  impulsive  component 
due  to  nearby  return  strokes  during  storm  activity  and  a  background  noise 
that  strongly  resembles  white  Gaussian  noise  (Ref  15).  The  Gaussian 
noise  is  due  to  the  large  number  of  storms  located  far  from  the  receiver 
plus  the  thermal  and  front  end  noise  of  the  receiver.  In  the  low 
frequency  channel,  the  thermal  and  front  end  noise  is  usually  negligi¬ 
ble  but  it  is  important  to  note  that  the  model  can  represent  either 
type  of  noise.  The  highly  impulsive  noise  component  is  assumed  to 
be  made  up  of  lightning  flashes  as  explained  in  Chapter  II.  Since 
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the  lightning  discharges  are  discrete  events  of  a  short  time  duration, 
it  is  reasonable  to  use  a  discrete  time  process  to  model  them.  Also, 
although  multiple  strokes  during  a  flash  may  be  correlated,  the  model 
will  be  designed  under  the  assumption  of  uncorrelated  strokes.  This 
is  done  to  preserve  simplicity  in  the  model  and  also  because  the 
correlation  will  have  little  effect  on  the  first  order  distribution  of 
the  noise.  The  assumption  of  uncorrelated  strokes  indicates  that  a 
Poisson  process  may  be  a  good  model  of  the  noise. 

Another  important  characteristics  of  lightning  that  will 
affect  the  model  is  that  the  ampliti.de  of  the  field  at  the  receiver 
is  different  for  each  stroke.  This  can  be  modeled  by  defining  a 
random  mark  for  each  Poisson  arrival  time.  The  random  amplitude 
will  be  modeled  as  a  zero  mean  Gaussian  random  variable. 

The  background  Gaussian  noise  can  be  modeled  by  a  similar 
process.  This  is  so  because  when  the  rate  of  the  Poisson  process  is 
high  relative  to  the  bandwidth  of  the  receiver,  the  process  observed 
at  the  output  of  the  receiver  will  be  indistinguishable  from  a 
Gaussian  process  (Ref  13). 

Therefore,  the  proposed  model  of  ARN  consists  of  two  separate 
marked  Poisson  processes.  One  process  has  a  high  rate  and  a  low 
variance  on  the  Gaussian  marks  to  model  the  white  Gaussian  background 
noise.  The  other  process  has  a  low  rate  and  a  relatively  high 
variance  on  the  Gaussian  marks  to  model  the  nearby  highly  impulsive 
ARN.  The  actual  variance  and  rates  of  the  two  processes  are  parameters 
to  be  calculated  for  each  specific  situation  to  be  modeled.  The 
details  of  the  parameters  choise  will  be  discussed  in  the  validation 
study. 
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The  resulting  noise  process  is  given  below. 


n(t)  =  l  a .6(t-t - )  +  l  bk6(t-uk) 

i=.o<>~  ~ 


(1) 


where 

n(t)  is  the  ARN  model 

are  the  independent  idnetically  distributed  Gaussian  marks 
on  the  background  process, 

t.j  are  the  Poisson  event  times  for  the  background  process, 

bk  are  the  independent  identically  distributed  Gaussian  marks 
on  the  impulsive  process, 

yk  are  the  Poisson  event  times  for  the  impulsive  process, 
and 

$(•)  is  the  Dirac  delta  function. 

In  the  next  section  section,  this  process  will  be  used  to  derive 
the  minimum  probability  of  error  processor  for  the  channel. 

Optimal  Receiver  Structure 

The  Bayes  optimal  processor  for  a  given  set  of  measurements  is 
the  processor  that  minimizes  the  expected  cost  associated  with  making  a 
decision  given  the  measurements.  In  this  report,  the  optimal  receiver 
will  be  a  Bayes  optimal  binary  hypothesis  receiver  with  a  symmetric  cost 
function  that  is  zero  for  a  correct  decision  and  some  positive  constant 
for  an  error.  This  minimum  probability  of  error  receiver  is  developed 
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in  many  texts  (Ref  17).  Given  a  set  of  measurements  in  the  form  of 
a  random  vector,  Eq.  (2)  expresses  the  decision  rule  that  will  mini¬ 
mize  the  probability  of  error. 


fR(R|Hj)  Hj  p 

W  <  *1  <2) 

~  H0 

where 

Pq  is  the  probability  of  hypothesis 

Pj  is  the  probability  of  hypothesis 

p 

;  is  the  random  measured  vector, 
and 

fR  (R | )  is  the  conditional  joint  density  of  the  measurement 
vector  given  that  signal  Si  is  sent  for  i = ( 0 , 1 ) . 
This  expression  is  called  the  likelihood  ratio  test.  The  joint 
densities  of  the  random  measurement  vector  indicate  the  necessity 
for  higher  order  statistics  of  the  channel  process  to  design  the 
optimal  receiver.  For  white  Gaussian  noise,  the  statistics  are 
known  from  the  first  order  density.  This  is  not  true  for  most  other 
processes.  Specifically,  in  the  case  of  atmospheric  radio  noise, 
the  joint  density  of  a  set  of  measurements  is  unsolved  for  the 
empirical  models  and  complex  for  the  proposed  poisson  process 
model . 


< 
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In  the  following  development,  a  minimum  probability  of 
error  receiver  for  the  low  frequency  channel  is  derived.  The  ARN 
model  will  be  used  to  represent  the  noise  in  the  channel. 

The  two  signal  hypotheses  that  the  receiver  must  decide 
between  are: 

H0:  r(t)  =  SQ(t)  +  n(t)  (3) 

H-| :  r(t)  =  S-j  (t)  +  n(t)  (4) 

where 

n(t)  is  the  ARN  model 

r(t)  is  the  measurement  waveform 

and 

Sg(t)  and  S^t)  are  marked  Poisson  processes  with  known  rate 

A  (t)  and  known  mark  density  f  (M). 

3  ~j 

In  the  following  discussion,  an  optimal  binery  receiver  using  the 
decision  rule  in  Eq.  (2)  will  be  derived. 

The  joint  density  of  the  marks,  event  times  and  the  number  of 
events  is  necessary  for  both  hypotheses  to  implement  the  optimal 
decision  rule.  The  joint  density  can  be  found  by  combining  the 
individual  Poisson  processes  into  a  single  Poisson  process.  This  is 
so  because  no  two  events  can  occur  simultaneously  in  a  Poisson  process 
(Ref  13).  The  resulting  density  of  a  mark  for  the  single  Poisson  pro¬ 
cess  will  be  a  mixture  density.  For  hypothesis  Hq,  the  mixture  den¬ 
sity  can  be  written  as. 
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I ti »Hq )  =  P[M^  came  from  the  background  process jt^jf  (M.) 
=  P[M^  came  from  the  impulsive  process  1 1^ 3f^ (M^ ) 

=  P[M.  came  from  the  signal  process |t.]f  (M.) 

1  1  !o  1 

It  is  implied  by  the  subscript  on  the  mark  that  the  mark  is 
associated  with  a  Poisson  event  time.  The  probability  that  the  mark 
came  from  a  particular  process  is  then  determined  by  the  relative  rates 
of  the  processes  at  the  event  time  associated  with  the  mark. 


'  VV+VV+xs.(ti>  ^ 


(5) 


The  marks  are  independent  of  each  other  so  the  joint  density  on  the  marks 
over  an  interval  of  T  given  that  N  evetns  occur  is  the  product  of  the 
individual  makr  densities.  The  event  time  statistics  for  a  Poisson 
process  is  given  by  Synder  (Ref  13). 


=  ft'N(t’N|HJ)fM(M|t,N,HJ) 


(6) 
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(7) 


T 

exp[-J  A_(t)+X.(t)+X  (t)dt] 

0  a  D  sj 

N  X,tti)f,(H1)nb(t1)f5(H1)HSj(t1)f!i(M1) 

1“  '  W+VV+xs.(V  J 

J 


By  writing  the  product  of  a  ratio  by  the  ratio  of  products  the  denomin¬ 
ator  will  cancel  and  yield 


I  X  ~  ^  J  *“J 


exp[/  X  (t)+X.(t)+X  (t)dt] 
0  a  b  Sj 


(8) 


The  resulting  likelihood  ratio  after  taking  the  natural  log  of  both 
sides  to  change  the  products  to  sums,  is  given  by. 


4 


( 
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(9) 


IP  T 

<  ln(/)  +  /  IX  ( t )— A  (t)]dt 

u  0  S1  s0 

Ho 

The  processor  therefore  turns  out  to  be  nonlinear.  This  means  that 
there  is  no  linear  correlation  receiver  that  can  implement  the 
processor. 

To  determine  if  this  processor  is  optimal  for  the  low  fre¬ 
quency  channel,  the  ARN  model  must  be  checked  for  its  representation 
of  the  noise.  In  the  next  chapter,  a  validation  study  will  confirm 
the  accuracy  of  the  ARN  model . 
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IV.  Verification  of  ARN  Model 


A  noise  model  must  be  proven  accurate  in  order  to  be  used 
with  confidence.  In  this  chapter,  the  first  order  statistics  of 
the  proposed  model  will  be  compared  to  those  of  measured  data. 

The  measured  data  is  in  the  form  of  the  amplitude  probability 
distribution  (APD)  of  the  envelope  of  ARN  detected  by  a  finite 
bandwidth  receiver.  APD  curves  measured  by  this  type  of  receiver 
are  presented  in  two  articles  that  will  be  addressed  in  this 
chapter.  They  are  (1)  K.  Furutsu  and  T.  Ishida  (Ref  3),  and 
(2)  CCIR  report  322  (Ref  22).  If,  through  a  choice  of  parameters 
for  the  model,  the  mean  square  error  (MSE)  is  small  between  the 
theoretical  APD  curve  given  by  Eq.  (39)  and  a  measured  APD  curve, 
then  the  model  will  be  assumed  to  be  a  good  representation  of  ARN. 

In  other  words,  if  the  first  order  statistics  of  the  model  are 
close  to  those  of  ARN,  then  the  higher  order  statistics  will  be 
assumed  to  be  close  also.  This  is  a  reasonable  assumption  because 
the  process  model  is  based  on  the  physics  of  the  noise  process 
as  was  discussed  earlier. 

This  chapter  consists  of  four  sections.  First,  (1)  the 
noise  model  will  be  used  as  an  input  to  a  mathematical  representation 
of  a  finite  bandwidth  envelope  detector.  The  resulting  APD  of  the 
envelope  at  the  output  of  the  receiver  will  be  used  for  validation. 
Next,  (2)  the  MSE  criterion  for  measuring  the  error  between  the 
measured  and  theoretical  APD  curves  will  be  discussed.  Then  (3) 
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the  Furutsu  and  Ishida  data  will  be  compared  to  the  theoretical 
APD  curves.  This  section  will  include  a  summary  of  the  Furutsu 
and  Ishida  article  plus  a  discussion  of  the  parameter  values  and 
their  relation  to  the  parameters  discussed  in  the  article.  Finally, 
(4)  the  data  will  be  compared  to  the  data  contained  in  CCIR  report 
322.  This  section  will  also  include  a  brief  discussion  of  para¬ 
meter  choice.  A  summary  of  CCIR  report  322  is  included  in 
Appendix  A  (Ref  6). 

In  the  next  section,  the  APD  of  the  output  of  a  quadrature 
receiver  is  derived  for  the  ARN  model. 

Verification  of  the  Point  Process 
Model 

To  verify  the  point  process  model,  it  will  be  used  as 
an  input  to  a  representation  of  the  receiver  used  in  the  CCIR 
report  322  and  the  K.  Furutsu  and  T.  Ishida  article.  Both 
receivers  are  envelope  detectors  with  a  finite  bandwidth.  A 
representation  of  this  type  of  receiver  is  shown  in  Figure  1. 

The  representation  will  be  referred  to  as  the  receiver  in  the 
following  developments.  The  receiver  consists  of  quadrature 
multipliers  followed  by  integrators  to  represent  the  bandpass 
filter  (Ref  13).  The  envelope  of  the  quadrature  components,  n^ 
and  n2»  is  then  detected  by  the  square  root  of  the  sum  of  squares 
of  the  quadrature  components.  In  this  section,  the  proposed  pro¬ 
cess  is  used  as  an  input  to  the  quadrature  receiver  described. 

The  resulting  statistics  of  the  envelope  at  the  output  will  then 
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be  compared  to  those  of  the  CCIR  report  322  and  the  Furutsu  article 
after  choosing  the  appropriate  parameters.  If,  through  a  choice  of 
parameters,  the  first  order  statistics  of  the  model  are  close  to 
those  of  the  measured  data,  the  higher  order  statistics  will  be 
assumed  to  be  close  also.  The  measure  of  closeness  or  fit  of 
the  data  will  be  discussed  later  in  this  section. 

The  input  process  to  the  receiver  is  defined  in  Eq.  (1). 

The  output  of  the  integrators  of  Figure  2  is  then, 


T 

n,  =  J  n(t)  cos  out  dt  (1°) 

-1  Q  ~  u 


T 

n 9  =  /  n(t)  sin  cont  dt  (H) 

0  ~  u 


Substituting  for  n(t)  in  Eq.  (16), 


2i 


T  «• 

=  J  l  a1-6(t-Ti )  cos  u^t  dt 


T  00 

+  /  l  b.  6(t-y.  )  cos  Wnt  dt 

0  k=-oo~K  ~K  u 


(12) 
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By  intercahnging  summation  and  integration  and  using  the  sifting 
property  of  the  delta  function, 


El 


l  ai  COS  0)nT.  +  l  b  cos  uvu. 
i=(r1  U~1  k=0~K  U~K 


(13) 


Where 


N  and  N,  are  the  number  of  events  which  occur  in  time  T 
~a  ~b 

for  the  a  and  b  processes  respectively.  Also, 


02 


l  a,-  sin  out.  +  l  b.  sin  u>np. 
i=0~i  u~i  k=0~K  u~k 


(14) 


The  problem  is  to  find  the  envelope  density  of  tne  n.^  and  n ^  quadrature 
components.  The  uk  and  x^  are  ordered  event  times  distributed  uni¬ 
formly  on  (0,T).  Usually,  T»1/ioq.  If  this  is  not  true  then  it 
must  be  assumed  that  uk  and  x^  are  uniformly  distributed  on 
(0,2nTr)  for  integer  n>0.  Also,  from  the  definition  of  the  marks, 
each  corresponding  pair  of  terms  in  n^  and  n 2  are  independent 
of  all  other  pairs  of  terms.  Defining  X  and  Y  to  represent  a 
typical  pair  of  terms  in  the  summations  for  n^  and  n^, 

X  =  a  cos  WqX  (15) 
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and 


Y  =  a  sin  u)qT. 

The  joint  density  of  X  and  Y  can  be  written  as, 
fXjY(X,Y|a)  =  ( R-  I  a  | )  (Ref  9:174) 

where 

R  =Vx2  +  Y2 


(16) 


(17) 


and 


6(R-|a|)  represents  the  probability  line  mass  fY  V(X,Y) 

A  5  I 

located  on  the  circle  X2+Y2=a2. 

The  characteristics  function  of  thejoint  density  on  X  and  Y  can 
now  be  written  as. 


V  (p»<J>|a) 
X.Y 


2n 


2tt  oo 


I  I  fy  Y(R’0la) 

0  0 


•  exp  [jRp  cos  (0-^)]RdRd9 


(18) 


where 
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6  =  tan-1  (y) 


This  is  the  joint  characteristic  function  for  the  random  variables 
X  and  Y,  even  though  this  integral  is  evaluated  in  polar  coordinates. 
By  the  sifting  property  of  the  delta  function. 


1  2lT 

(p,<j>|a)  =  j-J  exp[jap  cos  (6-0)]  de 
L  0 


=  JQ(ap) 


(19) 


where  JQ ( - )  is  the  zero  order  Bessel  function  of  the  first  kind. 
Each  corresponding  pair  of  terms  in  n-j  and  will  be  independent 
of  all  other  pairs  of  terms  in  the  summations.  This  allows  the 
joint  density  on  n^  and  n^  to  be  written  as  the  N-fold  convolu¬ 
tion  of  the  joint  density  on  each  pair.  In  the  characteristic 
function  domain,  this  convolution  can  be  written  as  a  product. 
Writing  the  joint  characteristic  function  of  n-j  and  as  the 
product  of  the  joint  characteristic  function  of  each  succeeeding 
pair  of  random  variables  in  the  summation  yields. 


l,'ni,n2(p’<,>lai’bk’Na’Nb) 


a  b 

n  Jn(a.p)  n  Jn(b  p) 
i=l  u  1  r=1  U  k 


(20) 


To  remove  the  conditioning  on  a.,  b.  ,  N  ,  and  N. ,  the  charac- 

*vK  ~  D 

teristics  function  can  be  averaged  over  the  statistics  of  a^,  b^, 
N  ,  and  respectively.  Since  each  term  in  the  first  and  second 
products  is  independent  of  all  other  terms,  the  averages  can  be 
done  for  each  term  individually. 


/  JQ(ap)  ffl(a)  da 

— oo 


00  1  2 

=  /  exp  da 

-00  V^TTO  2of 

a  a 

_p2o 2  p2o2 

=  exp  [~ IQ  (Ref  7:313]  (21) 


where  Ig  is  the  zero  order  modified  Bessel  function.  Performing  the 
same  average  over  the  statistics  on  eack  b^  term  leads  to  the  following 
result. 


2  2  2  2  H 
P  o  p  o  N 
a \ t  /  aii  a 


V?2(p’*lNa’V  =  [exp  (~T^)I0("4^)] 


-P2a?  p2a2  N. 
[exp(— 4-^)Iq(-^)]  b 


(22) 
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The  average  over  the  statistics  of  Nfl  and  Nfa  can  be  done  next. 


w  oo 


where 


,  (X  T)J'exp(-X  T) 

P^a=j)  =  — - Ji - §— 


(X.T)iexp(-A.  T) 

-  b  „  11 


T  is  the  integration  time  of  the  receiver. 
Performing  the  required  average  uields, 


WM)  ■  4* 


2  2  2  2 

P  °a  “P 
<  a\ _ / 


v  [XaTI0(~)exP(— 4~i)]Jexp(-XaT) 


22  2 
_  p  ”P  I 

•  exP[I0(— ^)exp( — ^ 


!bui|(AbT>'« 
n  i! 


exp(-AbT) 
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=  exp[XaT  I0(^i)exp(-^.)]exp(-XaT) 
p2o2  _  2a2 

•  exp[AbT  I0(-^)exp(-^)]exp(-AbT) 


f  PV  -pa2 

=  expj-T(Xa+xb)+T[XaI0(-^)exp(-^) 

P  2o2b  -p2a2  ) 

+  Vo(~4^)expH r^j 


(26) 


The  unconditional  joint  density  on  rij  and  n2  can  be  found  from  the 
unconditional  characteristic  function  through  the  transform  defined 
below. 


V.n.^’0)  =  TTTI  /  /  exP[~jRp  cos  (4,-0)] 
-1  ~2  (2tt)  0  0 


(27) 
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where 


R  is  now  defined  as 


nj+n 


7 

2 


and 


-1  n2 
0  =  tan  1(-^) 

nl 


,n2 


(R’e>  =  S  /  Jo(pR^, 


n  (p)pdp 

*  i*2 


(28) 


because  ^  (p,<}>)  is  independent  of  <b. 

~1~2 

Since  the  envelope  statistics  are  required,  it  is  necessary 
to  convert  the  joint  density  on  n^  and  to  a  joint  density  on  R 
and  6  where. 


R 


■i 


2  2 
~l+~2 


(29) 


and 


-1 

g  =  tan  J(^) 
1 
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fR  fl(R.e)  =  Rf  „  (R  cos  0,  R  sin  0) 

r'*v7  ni*''2 


=  S  4J0*pR)Vn2(p)pdc 


(30) 


The  cumulative  distribution  on  the  envelope  is  then  given  by, 


Rq  2-jj 

W  ‘  10  /0f?,e(R-e>dRde 


/  U  RJ0(pR)dR)p«.  „  (p)  dp 

o  o  u  5 


1  *~2 


Vi(pR0} 


Oi 


dp 


(Ref  7) 


{,R0Jl(pR0)*n1n2(P)  dp 


(31) 


The  amplitude  probability  distribution  (APD)  is  then  defined  by 

i~Fr(ro>  is» 
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A2 

p  °a, 


1  '  FR(R0}  =  1-/0RoJl(pR0)exP[-T(VXb)+T(XaIa(-r) 


2  2  2  2  2  2 
-p  a  p  aK  -p  a. 

exp(— — -)  +  AbI0(-^)exp(— ~)]dp 


(27) 


This  expression  can  be  numerically  integrated  for  values  of  Rg  to 
plot  the  APD  curve  for  a  set  of  given  parameter  values.  The  method 
of  integration  is  discussed  in  Appendix  B. 

This  APD  will  be  compared  to  measured  APD  curves  for  different 
parameter  values.  The  comparison  will  be  based  on  a  Mean  Square 
Error  (MSE)  distance  between  the  curves.  This  distance  is  discussed 
in  the  next  section! 

MSE  Criterion 

A  quantitative  measurement  of  error  between  APD  curves  must 
be  established  to  compare  the  theoretical  and  measured  APD  curves.  The 
APD  curve  is  a  plot  of  threshold  value  versus  the  probability  that  the 
envelope  of  the  noise  exceeds  the  threshold  value.  This  is  plotted 
on  Rayleigh  paper  which  forces  the  Rayleigh  distribution  to  plot  as 
a  straight  line.  The  criterion  to  be  used  to  measure  the  MSE  between 
curves  is  defined  below. 
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MSE 


(28) 


N 


2 


N 


Where 

MSE  is  the  sample  mean  square  error  between  the  APD  curves, 
d.j  is  the  linear  distance  between  the  curves  on  Rayleigh 

paper  for  a  given  threshold  value  as  given  by  Eq.  (29). 


di 


LOG 


rln<PMi\ 

10Lln(PTi)J 


(29) 


PMj  is  the  probability  that  a  measured  value  exceeds  the 
threshold  . 

Pyj  is  the  probabiliy  that  a  theoretical  value  exceeds  the 
threshold  . 
and 

N  is  the  number  of  samples  to  be  taken  from  the  curves. 

In  the  following  section,  the  Furutsu  and  Ishida  article  is 
discussed  and  two  measured  APD  curves  are  compared  to  theoretical 
APD  curves  of  the  ARN  model . 

Summary  of  K.  Furutsu  and  T.  Ishida  Article 

K.  Furutsu  and  T.  Ishida  have  developed  a  model  of  atmospheric 
noise  which  consists  of  a  Poison-Poison  process  (Ref  3).  A 
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Poison-Poison  process  is  a  Poisson  process  that  is  amplitude  modulated 
by  a  function  that  occurs  according  to  the  event  times  a  Poisson  pro¬ 
cess  of  another  rate.  They  based  this  model  on  the  natural  processes 
involved  in  lightning,  the  major  source  of  atmospheric  noise.  The  pre¬ 
discharge  (stepped  leader)  consists  of  many  random  impulses  of  the 
order  of  duration  of  1  usee,  each  of  which  occur  successively  in  the 
time  interval  of  25  to  100  usee.  This  indicates  that  the  pre¬ 
discharge,  which  would  occur  in  time  according  to  a  Poisson  distribu¬ 
tion,  forms  a  wave  packet  of  Poisson  noise.  In  the  paper,  these 
packets  have  an  envelope  of  a  growing  exponential  which  reaches  a 
finite  value  and  then  returns  to  zero. 

The  process  is  then  validated  by  comparing  the  theoretical 
first  order  density  to  the  measured  first  order  density  of  the  output 
of  an  envelope  detector  with  a  finite  bandwidth.  The  receiver  detects 
the  amplitude  waveform  with  a  bandpass  filter  and  then  measures  the 
amplitude  distribution  of  the  envelope. 

Two  cases  are  used  for  the  validation  study.  One  case  assumes 
that  the  noise  is  created  by  a  strong  local  source  at  a  short  distance 
and  that  other  sources  are  continuously  distributed  far  distant  so 
that  the  noise  strength  is  much  lower  for  the  distant  sources  starting 
at  some  distance  from  the  receiver  and  extending  to  infinity.  The 
amplitude  of  the  noise  varies  inversely  with  the  distance  from  the 
receiver.  It  is  clear  from  their  study  that  there  is  strong  correla¬ 
tion  between  their  model  and  the  measured  ARN. 
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Comparison  to  Furutsu 
and  Ishida  Data 

To  begin  the  comparison,  the  parameters  of  the  theoretical 
model  developed  in  Chapter  III  must  be  chosen.  Since  the  model  that 
Furutsu  and  Ishida  developed  is  related  to  the  Poisson  process,  the 
parameters  that  they  specify  for  the  measured  APD  curves  are  directly 
applicable  to  the  model  of  Chapter  III.  The  parameters  of  the  two  APD 
curves  in  Figs.  2  and  3  are  given  in  Table  I.  The  only  parameter  that 
had  to  be  experimentally  determined  was  the  rate  of  the  background 
process.  It  was  found  that  a  rate  of  about  1.5  to  1.8  times  the  band¬ 
width  of  the  filter  was  sufficient  to  produce  good  agreement.  In  the 
Furutsu  and  Ishida  curves  the  envelope  threshold  (plotted  vertically 
for  the  two  plots)  was  in  terms  of  a  ratio  of  the  envelope  threshold 
to  the  standard  deviation  of  the  background  noise.  Therefore,  the 
background  noise  for  the  point  process  model  is  chosen  to  be  1.  This 
allows  both  curves  to  be  easily  plotted  on  the  same  scale.  The  other 
parameters  are  defined  in  terms  of  two  other  ratios.  The  ratios  are 
(1)  the  rate  of  the  impulsive  process  divided  by  the  bandwidth  of  the 
IF  filter  and  (2)  the  standard  deviation  of  the  impulsive  noise 
divided  by  the  standard  deviation  of  the  background  noise.  The  second 
ratio  is  expressed  in  terms  of  decibels  where, 

value  in  Db  =  20  log10  (ratio) 

In  addition  to  these  parameters,  the  integration  time  was  chosen  by 
assuming  that  the  integrator  is  approximately  an  ideal  bandpass  filter 
and  that  the  double  sided  bandwidth  is  2/T.  T  is  the  integration 
time.  Actually,  the  integrators  have  a  sine  shaped  respone  in  the 
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Fig.  3.  Comparison  of  Measured  and  ARN  Model 
APD  Curves  (MSE=.0l6) 
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frequency  domain  but  close  agreement  was  found  using  this  approxima¬ 
tion. 


TABLE  I 


PARAMETERS  FOR  THE  ARN  MODEL 


Figure 

Xa 

Xb 

°a 

ab 

T 

3 

2500.0 

30.4 

1.0 

158.5 

.00132 

4 

2000.0 

72.0 

1.0 

63.0 

.00167 

5 

350.0 

20.0 

.11 

2.0 

.01 

It  is  found  that  the  theoretical  curves  from  the  point  process 
developed  in  Chapter  III  are  a  very  close  fit  to  those  of  the  measured 
data  in  the  Furutsu  and  Ishida  report.  This  is  significant  because  the 
data  was  taken  in  the  daytime  with  nearby  storms.  This  indicates  that 
the  known  correlation  between  strokes  for  nearby  storms  is  adequately 
represented  in  the  first  order  statistics  derived  from  the  model. 

Next,  a  similar  comparison  to  the  CCIR  report  322  data  will  be  made. 

Comparison  to  CCIR  Report  322  Data 

In  CCIR  Report  322,  measured  ARN  data  from  16  receiving  sta¬ 
tions  located  around  the  world  is  collected  and  averaged.  A  summary 
of  the  report  is  given  by  Hettinger  (Ref  6)  in  Appendix  A.  Several 
APD  curves  are  given  in  the  report  and  each  is  specified  by  only  a 
Vj  ratio  which  is  a  measure  of  the  impulsiveness  of  the  noise.  The 
definition  of  this  ratio  is  given  in  Eq.  (A4).  To  relate  the  param¬ 
eters  of  the  point  process  model  to  this  ratio  would  require 
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computations  of  ensemble  averages.  The  averages  were  not  computed  due 
to  the  time  constraints  of  this  project  (they  require  numerical 
integrations).  Since  the  parameters  of  the  model  could  not  be 
directly  related  to  the  Vd  ratio,  each  parameter  was  determined  by 
the  parameter's  effect  on  the  resulting  theoretical  APD  curve.  The 
rate  of  the  background  proces  was  set  to  1.75  times  the  bandwidth, 
which  was  determined  in  the  Furutsu  and  Ishida  comparison.  Also,  an 
article  by  Spauling  et  al.  (Ref  15)  indicated  that  the  lower  portion 
of  the  curve  is  determined  by  the  background  noise  and  the  upper  por¬ 
tion  is  determine  by  the  impulsive  noise.  Therfore,  a  change  in  the 
variance  of  the  background  process  marks  produced  a  fit  to  the  lower 
portion  of  the  curve.  The  upper  portion  of  the  curve  was  then  matched 
by  a  variation  of  both  the  rate  and  variance  of  the  impulsive  process. 
Although  this  may  seem  arbitrary,  for  given  ARN  the  rate  of  the 
impulsive  process  can  be  determined  by  observing  the  number  of  large 
fluctuations  in  the  envelope  of  the  noise  in  a  given  time.  The 
variance  can  then  be  used  to  adjust  the  power  of  the  model  so  that  it 
matches  the  noise  power. 

The  plot  shown  in  Figure  4  shows  that  for  the  chosen  param¬ 
eter  values  the  model  acurately  represents  the  first  order  statistics 
of  ARN  as  given  for  Vd=8  in  the  CCIR  report. 

The  point  process  model  is  therefore  found  to  be  versatile  in 
its  ability  to  represent  the  first  order  statistics  of  ARN.  The  para¬ 
meters  of  the  model  can  be  chosen  to  represent  a  variety  of  situa¬ 
tions,  of  which  three  are  shown  here. 

It  is  interesting  to  note  that  if  other  types  of  noise  that 
are  not  correlated  with  the  ARN  or  each  other  are  significant,  the 
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model  allows  as  many  marked  point  processes  as  are  necessary  to  be 
added  wthout  increasing  the  complexity  of  the  analysis  done  in  this 
report.  This  is  another  useful  feature  of  the  model.  In  the  next 
chapter,  some  suggestions  will  be  made  for  further  studies  along  with 
some  conclusions  about  the  point  process  model. 
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V.  Conclusions  and  R commendations 


In  conclusion,  a  sum  of  marked  Poisson  processes  is  a  good 
representation  of  ARN.  The  ARN  model  developed  in  this  report  was  the 
sum  of  two  marked  point  processes,  each  with  zero  mean  Gaussian  marks. 
Even  with  this  simplified  structure,  the  model  was  found  to  be  in  good 
agreement  with  the  measured  first  order  distribution  of  ARN.  The 
simplified  structure  constitutes  a  tractable  model  that  can  be  traced 
through  a  quadrature  receiver,  and  results  in  the  design  of  a  Bayes 
optimal  binary  receiver  for  the  channel. 

It  is  recommended  that  further  studies  of  this  model  be 
directed  to  two  areas.  These  are  (1)  a  more  thorough  validation 
should  be  done  with  raw  data.  This  validation  should  involve  higher 
order  statistics  of  the  process  and  should  determine  the  useful  range 
of  the  model.  (2)  Performance  evaluations  should  be  done  for  both 
optimal  receivers  designed  for  the  low  frequency  channel  and  the 
suboptimal  receivers  presently  being  used  in  the  low  frequency  channel. 
These  performance  evaluations  should  be  done  using  the  marked  Poisson 
process  model  as  the  input  noise  to  the  receivers.  This  will  allow 
bounds  on  the  performance  of  the  receivers  to  be  determined  along  with 
a  comparison  of  receivers  presently  in  use. 
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Appendix  A:  Summary  of  CCIR  Report  322 
by  Hettinger  ( Ref  67 


A  large  quantity  of  ARN  data  has  been  gathered  by  the  National 
Bureau  of  Standards.  This  data  includes  worldwide  time  averaged  noise 
power  measurements,  and  amplitude  probability  distribution  curves  for 
the  time  varying  noise  envelope.  This  work  is  presented  in  condensed 
form  in  CCIR  Report  322  (Ref  22).  The  block  diagram  of  the  receiver 
used  in  these  measurements  is  presented  in  Figure  Al. 

The  worldwide  measurements  of  ARN  power  are  derived  from  mea¬ 
surements  made  by  sixteen  receiving  stations  throughout  the  world. 
Figure  A2  shows  the  location  of  these  stations.  The  average  noise 
power  available  at  the  antenna  is  relatively  constant  for  observa¬ 
tion  intervals  up  to  four  hours.  The  measurements  made  during  these 
four-hour  periods  are  relatively  constant  from  day-to-day  for  periods 
up  to  three  months.  Therefore,  the  noise  power  measurements  in  CCIR 
Report  322  are  plotted  for  the  four  seasons  of  the  year,  and  each 
seasonal  section  contains  six  plots  of  the  four  hour  time  blocks  in 
a  day  (Ref  22).  Figure  A3  shows  a  typical  average  noise  power  plot 
from  CCIR  Report  322.  The  average  noise  figure,  F  ,  shown  in  Figure 
A3  is  computed  from  Eq  (Al). 


Ffl  =  10  Log 
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Fig.  A3.  A  Typical  Plot  of  Worldwide  Average 
Noise  Power  as  in  CCIR  Report  322 
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where 


F  =  the  average  noise  figure  in  db 

d 

K  =  Boltamn's  constant 

Tq  =  Standard  temperature  of  280°  Kelvin 
E(t)  =  The  time  varying  envelope  waveform 
T  =  The  observation  time  in  sec 
B  =  the  receiver  bandwidth  in  Hertz 

All  of  the  noise  power  measurements  in  CCIR  Report  322  are 
made  by  receivers  with  a  200  Hz  effective  bandwidth.  Figure  A4 
shows  an  example  of  an  empirically  derived  plot  for  converting  the 
noise  power  measured  at  1  MHz  to  the  average  noise  powerat  a  dif¬ 
ferent  frequency.  A  different  noise  power  conversion  plot  is  required 
for  each  worldwide  plot  of  noise  power. 

CCIR  Report  322  also  contains  a  plot  of  the  amplitude-proba¬ 
bility  distribution  (APD)  of  the  ARN.  The  APD  curve  is  a  plot  showing 
the  percenage  of  the  time  that  the  time  varying  envelope  voltage, 
at  the  output  of  an  envelope  detector,  exceeds  some  threshold.  The 
curve  is  plotted  parametrically  with  V^,  which  is  defined  by  Eqs  (A2), 
(A3),  and  (A4).  The  Vd  ratio  is  a  rough  measure  of  the  impulsiveness 
of  the  noise;  the  more  impulsive,  the  larger  the  ratio. 
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Fig.  A4.  Typical  Noise  Power  Conversion  Plot 
as  in  CCIR  Report  322 
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Eav  -  [{  J  «t)dt] 
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The  shapes  of  the  APD  curves  depend  only  on  the  ratio, 

rather  than  the  absolute  value  of  E_  or  E..f  (Ref  1:1).  Thus,  all 

of  the  APD  curves  with  a  given  ratio  can  be  collapsed  into  a  sngle 

curve  by  normalizing  the  incoming  waveform  by  Erms>  such  that  the 

resulting  waveforms  all  have  power  equal  to  1  watt  (Ref  22:8-9). 

Since  the  waveforms  have  been  normalized  by  E  .  the  envelope  vol- 

rms 

tage  threshold  E-j.  must  also  be  normalized.  Therefore  the  ordinate 
of  the  APD  plots  is  presented  in  terms  of  a  threshold,  A,  which  is 
defined  in  Eq  (15) . 


A  =  20  L0G1q 


(A5) 


where 

Ej  =  envelope  voltage  threshold. 

The  APD  plot  shown  in  Figure  A5  represents  the  amount  of  time  the 
envelope  waveform  exceeds  a  threshold  for  a  given  V^. 
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Both  the  noise  figure,  F  ,  and  the  APD  curves  are  dependent 

•  a 

on  the  receiver  bandwidth  (Ref  15).  The  receiver  used  by  the  National 
Bureau  of  Standards  has  an  effective  noise  power  bandwidth  of  200  Hz. 

The  APD  curves  in  Figure  A5  can  be  modified  to  apply  to  receivers  of 
bandwidth  other  than  200  Hz.  This  can  be  done  by  selecting  a  curve  with 
a  different  Vd.  Figures  A6  and  A7  can  be  used  to  convert  a  Vd  ratio 
measured  by  a  receiver  of  one  bandwidth  to  that  of  a  different  band¬ 
width  (Ref  15). 

The  APD  curves  in  CCIR  Report  322  are  composed  from  measurements 
made  in  different  seasons  of  the  year  and  times  of  the  day.  Due  to 
the  long-term  variations  that  occur  with  time,  the  curves  can  only  be 
considered  an  approximation  to  an  APD  curve  that  might  be  measured  over 
any  short  period  of  time.  Crichlow  does  not  believe  this  introduces 
an  appreciable  error  in  computations  using  the  APD  curves  (Ref  1:10). 

All  of  the  measurements  made  in  CCIR  Report  322  relate  to  the 
statistics  of  the  time  varying  envelope  waveform  of  the  noise.  If 
the  envelope  noise  process  is  modeled  as  a  stochastic  process  that  is 
ergodic  in  both  the  mean,  and  autocorrelation,  then  the  measure¬ 
ments  of  average  noise  power  and  the  average  envelope  amplitude  from  CCIR 
Report  322  are  representative  of  the  second  and  first  moments  of 
the  first  order  probability  density  function  of  the  nofse  process. 

If  the  process  is  ergodic  in  the  mean,  then  it  is  also  ergodic  in  dis¬ 
tribution  (Ref  9:328-332).  Therefore,  the  measured  APD  curves  can 
be  related  to  the  first  order  cumulative  distribution  function  of  the 
envelope  noise  process.  The  assumption  of  ergodicity  over  at  least 
four-hour  periods  of  time  is  required  to  make  the  noise  measurements 
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Bandwidth  Ratios  B2/B 


Fig  A7.  APD  Bandwidth  Conversion  Plot  2 


from  CCIR  Report  322  a  meaningful  measure  of  the  statistical  properties 
of  the  random  process. 


Appendix  II :  Numerical  Integration  of 
Equation  38 

The  purpose  of  this  appendix  is  to  discuss  the  numerical 
integration  of  FR(RQ)  given  at  the  end  of  Chapter  IV.  This  equa¬ 
tion  is 


fr(r0)  =  exp  L-T( Afl+Ab) 

Ro 

•  R0  {/l(pR0)exp[T(XaI0(Za)exp('Za) 

*  WZb,exp(‘Zt>>J  dp  (B.l ) 
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The  exponential  form  of  IQ( Z)  given  in  Eq.  (B.2)  shows  that 
for  large  p,  the  exponent  will  approach  zero. 


Vz>  - 


v^ttZ 


-(1 


(Ref  17) 


(B.2) 


for  Z»1 . 

It  then  follows  that  the  exponential  term  approaches  the 
value  1  and  eventually  has  little  effect  on  the  integral  as  p  gets 
large.  Therefore,  the  integral  is  done  in  two  sections. 

The  first  section  of  theintegral  was  done  using  the  tra¬ 
pezoidal  rule.  This  section  is  ended  when  the  exponent  becomes 
less  than  a  set  threshold  value.  Since  for  small  X,  exp(X)=l+X 
so  the  threshold  is  directly  related  to  the  error  in  the  approxi¬ 
mation  that  exp(X)=l.  The  threshold  value  used  for  the  given 
curves  was  .0001  which  allowed  the  integral  to  be  evaluated  in  a 
reasonable  amount  of  time.  The  bessel  functions  were  evaluated 
using  a  routein  supplied  by  the  computer. 

The  final  section  is  then 


/  RQJi(pR0)dp 
P1 


where  pj,  is  the  last  value  of  the  first  integration.  By  letting 
u=pRq,  dp=du/Rg  so  that  the  integral  is  reduced  to 
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L  /R  Jl(U)dU  =  J0(^} 

P^Rq  0 


Therefore,  the  integral  can  be  defined  as. 


W  =  exP[-T(xa+?'b)3[Ro5;i+J0(R^)] 


(B.3) 


where  is  the  integral  to  . 
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